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White-rot and brown-rot flingi employ different mechanisms to degrade lignocellulose. These fungi are not monophyletic 
and even alternate in their common lineage. To explore the reason for this, seventy-six ligninolytic peroxidases (LPs), 
including 14 sequences newly identified from available basidiomycetous whole-genome and EST databases in this study, 
were utilized for phylogenetic and selective pressure analyses. We demonstrate that LPs were subjected to the mixed 
process of concerted and birth-and-death evolution. After the duplication events of original LPs, various LP types may 
originate from mutation events of several key residues driven by positive selection, which may change LP types and even 
rot types in a small fraction of wood-decaying fiingi. Our findings provide preliminary insights into the cause for the 
alternation of the two fungal rot types within the same lineage. 
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Introduction 

Wood-decaying flingi could release nutrients, particularly 
carbon, from wood, and thus play an essential role in the 
function of forest ecosystems. For this reason, they attracted 
widespread attention of forest ecologists (Lonsdale et al. 
2008) and other scientists representing different disciplines, 
such as fiingal taxonomy (Lutzoni et al. 2004), forest 
pathology (Asiegbu et al. 2005), and biotechnological 
applications (Cohen et al. 2002). Two main rot types of 
wood are known, namely, white-rot and brown-rot (Worrall 
et al. 1997). The former is caused by white-rot fiingi, which 
can almost totally decay lignocelluloses, including Hgnin, 
parts of cellulose, and hemicellulose. The latter decay is 
caused by brown-rot fiingi, which only selectively remove 
cellulose and hemicellulose (Ryvarden 1991). White-rot 
and brown-rot fiingi employ two completely different 
mechanisms to fiilly and partially degrade components of 
lignocellulose, and the most important difference between 
them is whether they have capacity to degrade hgnin, one of 
the most abundant natural biopolymers. Therefore, they 
might have different evolutionary origins. Actually, the 
degradation capacity has been used as an important taxo- 
nomic character in Basidiomycota for more than 20 years 
(Redhead and Ginns 1985). Surprisingly, in the phylogeny 
of Basidiomycota, inferred from the combination of nuclear 
and mitochondrial small-subunit ribosomal DNA data, 
white-rot and brown-rot fiingi are not separated from each 
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Other (Hibbett and Donoghue 2001). Moreover, brown-rot 
fiingi have six independent origins within white-rot fimgal 
branches and inversely the white-rot fungus, Grifola fron- 
dosa (Dicks.) Gray, is situated on a branch of brown-rot 
fiingi (Hibbett and Donoghue 2001). Based on ancestral 
state reconstruction, Hibbett and Donoghue (2001) pro- 
posed that the cause for the bidirectional transformation 
between white-rot and brown-rot fungi is that the genes 
encoding lignin-degrading enzymes of white-rot fungi are 
also retained in the genomes of brown-rot fungi, and the 
difference is caused just by the change of gene expression 
patterns. This explanation, however, does not take the 
available data on expression systems into consideration. 
Moreover, the genome of Postia placenta (Fr.) M.J. 
Larsen and Lombard, a brown-rot fungus, lacks genes 
encoding true lignin-degrading enzymes (Martinez et al. 
2009). Therefore, it is still not fitUy understood that how 
white-rot and brown-rot fimgi alternate in their lineage from 
the view of lignin-degrading enzymes, of which two kinds 
are known, namely, ligninolytic peroxidases (LPs) and 
laccases. In this study, we focused on the LPs, because 
their reaction mechanism has been intensively studied 
(Hammel and Cullen 2008). 

LPs are extracellular fungal peroxidases, belonging to 
Class-II of the plant/fiingal/bacterial peroxidases superfam- 
ily, which also includes two other classes: Class-I, intracel- 
lular peroxidases and Class-Ill, extracellular plant 



© 2014 The Author(s). Published by Taylor & Francis. 

This is an open-access article distributed under the tenns of the Creative Commons Attribution License http://creativecommons.Org/licenses/by/3.0/, which permits unrestricted 
use, distribution, and reproduction in any medium, provided the original work is properly cited. The moral rights of the named author(s) have been asserted. 



30 L.-W. Zhou et al. 



peroxidases (Conesa et al. 2002; Passardi, Theiler, et al. 

2007) . All members of the three classes share the same 
heme moiety and are believed to be evolutionary closely 
related (Passardi, Bakalovic, et al. 2007; Morgenstem et al. 

2008) . Phylogenetic analyses of this superfamily are mainly 
focused on Class-I (Zamocky 2004; Passardi, Zamocky, 
et al. 2007), even if Passardi, Bakalovic, et al. (2007) 
explored the whole relationship of this superfamily. More 
recently, Morgenstem et al. (2008) focused on Class-II, but 
concentrated mainly on the description of LP diversity, and 
whether the various types of LPs were monophyletic or 
polyphyletic rather than on the molecular dynamics diversi- 
fying the types. Until now, there is no comprehensive and 
in-depth phylogenetic analysis of LPs, which are distin- 
guished from other members of this superfamily by the 
capacity to degrade lignin. Three well-known LP types 
exist in white-rot fungi: hgnin peroxidases (LiP), manganese 
peroxidases (MnP), and versatile peroxidases (VP). Various 
white-rot fungi can produce either one, two, or all of them. 
LiP was first described in the growth medium of 
Phanewchaete chrysosporium Burds (Kirk and Farrell 
1987). It was characterized by the capacity to oxidize high 
redox potential aromatic substrates, such as non-phenolic 
lignin structures, for its tryptophan residue (W171 referred 
to as LiPA of P. chiysosporium; Doyle et al. 1998). 
Meanwhile, MnP was identified in the same species (Gold 
et al. 1989). In comparison to LiP, MnP is unable to oxidize 
high redox potential substrates for lack of the key tryptophan 
residue (Hammel and CuUen 2008). However, MnP can 
oxidize the relatively low redox potential phenolic lignin 
structures by the oxidation of Mn^^ by a manganese-binding 
site (E35xxxE39...D179 referred to as MnPl of 
P. chrysosporium; Wariishi et al. 1991; Sundaramoorthy 
et al. 1994). Because VP is absent in the best studied ligni- 
nolytic fungus P. chrysosporium, it was not discovered until 
the late 1990s, in Pleurotus eryngii (DC.) Quel. (Camarero 
et al. 1999). VP is known for sharing the oxidation sites of 
both LiP and MnP. In addition to the three types, a fourth 
type of LPs, which has lost all or part of the key residues 
essential in the oxidization reactions, was first found in 
Coprinopsis cinerea (Schaeff.) Redhead et al, by 
Baunsgaard et al. (1993). Although the tertiary structure is 
similar to that of the other types (Larrondo et al. 2005), it still 
needs to be investigated if it can efficiently degrade lignin for 
the lack of the key residues. We called it 'CII' short for 'other 
class II peroxidase' in this study following the nomenclature 
of PeroxiBase (Passardi, Theiler, et al. 2007). LPs often have 
several isozymes, such as the 10 LiPs and 5 MnPs in P. 
chrysosporium (Martinez et al. 2004; Vanden 
Wymelenberg et al. 2006). This may reflect the specificity 
of the LPs in response to the diversification of lignin 
(Passardi, Bakalovic, et al. 2007). 

By the comparison of whole genome sequences of the 
phylogenetically closely related white-rot fungus P. chry- 
sosporium and brown-rot fimgus Postia placenta. 



Martinez et al. (2009) suggested that brown-rot fimgi 
were shift from white-rot fungi by losing the capacity to 
degrade lignin, particularly by losing LPs. However, the 
phylogeny was constructed using only a few LPs, and 
their origin in white-rot fiingi was not shown. It is pro- 
posed that the LPs originated from a cytochrome c perox- 
idase (CcP), belonging to the Class-I of the plant 
peroxidases superfamily, which changed and became 
able to degrade the novel appearing high redox potential 
polymer lignin (Passardi, Bakalovic, et al. 2007). 
However, the reason that different fimgal species express 
different types of LPs and the evolutionary order of LPs 
were still unknown. Recently, Ruiz-Duenas et al. (2013) 
explored the diversity of LPs from 10 sequenced genomes, 
but all resulted from Polyporales. Therefore, these preli- 
minary studies could not convincingly explain the evolu- 
tionary origin and the shift of white-rot and brown-rot 
fungi within the basidiomycetous lineage. 

By searching the available basidiomycetous whole-gen- 
ome and EST databases, we extended the annotated data 
set from PeroxiBase (Passardi, Theiler, et al. 2007) and 
assembled the most complete and nonredundant LPs with 
the information of gene structures. We explored the phylo- 
geny of LPs and elucidated their evolutionary dynamics in 
both white-rot and brown-rot fungi taking into considera- 
tion both the types of LPs and fungal rot. The results 
presented below can provide prehminary insights into the 
nature of the alternation of white-rot and brown-rot fungi 
in their fimgal lineage at the molecular level. 



Materials and methods 
Data mining 

All sequence data used in this study were retrieved from 
PeroxiBase (Passardi, Theiler, et al. 2007), JGI (the DOE 
Joint Genome Institute website: http://www.jgi.doe.gov/) 
and NCBI (National Center for Biotechnology 
Information: http://www.ncbi.nlm.nih.gov/). LP amino 
acid sequences were obtained ixom PeroxiBase, the latest 
update on 22 September 2011 for LPs, and partial ones 
were completed using the corresponding genome 
sequences in JGI. To avoid redundancy, if there are several 
same LPs from various strains of one fimgal species, only 
one LP amino acid sequence was selected for further 
analyses, and all incompletely processed transcripts 
(Macarena et al. 2005) were discarded. In addition, certain 
amino acid sequences were edited manually according to 
their genome sequences. Ip gene sequences (including 
both exons and introns) were retrieved from the JGI or 
NCBI databases, if available. Six CcPs were selected as 
outgroups for the phylogenetic analysis. 

To identify novel LP sequences, all available basidio- 
mycetous genome sequences of both white-rot and brown- 
rot fungi in JGI and NCBI were screened using tBlastN 
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with default settings and different LP amino acid 
sequences as queries. All queries with high sequence 
similarity produced similar hits. Exon-intron boundaries 
of the identified genomic sequences were recognized 
based on the alignments of ESTs and gDNA and the 
conserved intron splicing rule. The obtained coding 
sequences were translated into amino acid sequences 
using the ExPASy (Expert Protein Analysis System) trans- 
lation tool (http://www.expasy.org/tools/dna.html). The 
conserved residues mainly in charge of redox fianctions 
were detected through sequence alignments using 
ClustalX 1.81 (Thompson et al. 1997) with default settings 
and comparisons of the tertiary structure using 3D- 
JIGSAW (Contreras-Moreira and Bates 2002) with pre- 
viously verified LPs (reviewed by Martinez (2002)). 
Retrieved sequences without conserved heme-binding 
sites were discarded; however, those lacking key residues 
for oxidizing substrates, but with conserved tertiary struc- 
tures, were included to make the analysis complete. 

Sequence alignments and phylogenetic analyses 
SignalP 3.0 online server (Bendtsen et al. 2004) was used 
to predict the signal peptides (SPs) of all LPs and six CcPs 
using the neural network model with default settings. Two 
protein data sets, one is all LPs and six CcPs, the other is 
those without SPs, were used in subsequent phylogenetic 
analysis. The two data sets were aligned with ClustalX 
1.81 using default settings (Gonnet series weight matrix, 
gap opening: 10, and gap extention: 0.2) and manually 
edited, respectively. 

Three methods were performed for phylogenetic ana- 
lyses of the two data sets. MEGA 3.1 (Kumar et al. 2004) 
was used to construct neighbor-joining (NJ) trees based on 
pair-wise deletion of gaps/missing data and p-distance 
matrix of amino acids model with 1000 bootstrap repli- 
cates. The maximum likelihood (ML) trees with 100 boot- 
strap replicates were constructed using PhyML 3.0 
(Guindon and Gascuel 2003) with the best-fit evolutionary 
model (WAG + I + G + F) selected by ProtTest 2.4 
(Abascal et al. 2005) according to the AIC criterion. The 
maximum parsimony (MP) trees with 1000 bootstrap 
replicates were constructed using PAUP* 4.0b 10 
(Swofford 1998). All characters have equal weight, and 
gaps are treated as missing. 

Gene structure analysis 

For all known Ip introns located in neither 5' nor 3' 
untranscripted region (UTR), only the gene sequences 
located between the translation start and stop codons 
were studied. All 66 full-length amino acid sequences 
with detailed gene structures were aligned with ClustalX 
1.81 using default settings. Position and phase of introns 
were recognized based on the comparison of boundary 



exons and their corresponding amino acid sequences, 
and mapped onto the matrix of multiple alignments. The 
intron phase was set as 0, if an intron intervenes between 
two consecutive codons, 1 if an intron intervenes between 
the first and the second codon nucleotides, and 2 if an 
intron intervenes between the second and the third codon 
nucleotides, as described in Sanchez et al. (2003) and 
Zhou et al. (2008). To visualize the comparison of gene 
structures clearly, the introns were renamed using capital 
letters in alphabetical order corresponding to their appear- 
ing order in the above matrix, and the LPs fi'om the same 
groups were grouped together manually. 



Selective pressure estimation 

The CODEML program in PAML 4.2 (Yang 2007) was 
used to estimate the selective pressure on LPs. Following 
the recommendation of Anisimova et al. (2002), three 
pairs of site models, namely, MO (one-ratio) vs. M3 (dis- 
crete), Mia (NearlyNeutral) vs. M2a (PositiveSelection), 
and M7 (beta) vs. M8 (beta and co) (Yang et al. 2000; 
Wong et al. 2004), were employed to detect positive 
selection. The values of nonsynonymous/synonymous 
substitution rate ratio ((o ~ dN/dS) < 1, = 1, and > 1 
mean purifying selection, neutral evolution, and positive 
selection, respectively. Therefore, only when the co value 
of model M3, M2a, or M8 was higher than one, a like- 
lihood ratio test (LRT; Nielsen and Yang 1998) was used 
to test the results of the above three nested models. The 
sites under positive selection were also identified by both 
naive empirical Bayes (NEB) and Bayes empirical 
Bayes (BEB). 



Results 

Identification and distribution of LPs 

A total of 76 LP amino acid fiiU-length sequences (in 
Supplementary File 1) were identified, including 14 
novel ones, of which 4 and 10 were completed and pre- 
dicted from genome sequences, respectively. The corre- 
sponding genomic data of 66 LPs were successfiilly 
retrieved. Relevant infonnation is listed in Table 1. The 
76 LPs fell in three fungal orders: Russulales, Agaricales, 
and Polyporales, as described by Passardi, Bakalovic, 
et al. (2007). Of these LPs, 20 are LiPs exclusively 
found in the Polyporales, 11 are VPs found in both the 
Agaricales and Polyporales, 37 are MnPs, and 8 are CIIs, 
both of which are found in all three orders. The alignment 
of 24 selected LPs, including CII, MnP, LiP, and VP, was 
presented to show the key residues essential in the oxidi- 
zation reactions (Figure 1). 

All other available whole-genome sequences in JGI 
and NCBI Irom 11 species/variant species, which repre- 
sent six basidiomycetous orders (Table 2), were also 
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Figure 1. Partial sequence alignment of 24 selected LPs. Their code names are tiie same as those in Table 1. Dot lines mean the omitted 
amino acid residues. The numbers in front of and behind the residues, respectively, refer to the positions of the first and the last residues 
in its fijU-length amino acid sequence. The key residues, ExxxE...D and W, essential in the oxidization reactions are shaded gray. The 
LPs with ExxxE...D, W, both and neither are MnP, LiP, VP, and CIL respectively. 



Table 2. Species without LPs and their taxonomic position. 


Species 


Order 


Species 


Order 


Ciyptococcus bacillisporus 
Cryptococcus neoformans var. grubii 
Cryptococcus neoformans var. neoformans 
Malassezia globosa 
Malassezia restricta 
Melampsora laricis-popiiUna 


Tremellales 

Tremellales 

Tremellales 

Malasseziales 

Malasseziales 

Pucciniales 


Moniliophthora perniciosa 
Puccinia graminis f.sp. tritici 
Schizophyllum commune 
Sporobolomyces roseus 
Ustilago maydis 


Agaricales 

Pucciniales 

Agaricales 

Sporidiobolales 

Ustilaginales 



screened using tBlastN, but no positive hits according to 
our filtration criterion (see Data mining for details) were 
retrieved, indicating that these species perhaps harbor 
no LPs. 

Phylogenetic relationship of LPs 

In the alignment of data set with SPs, 507 characters with 
376 infomiative ones exist. NJ, ML, and MP trees 
(Supplementary Figure SI) were constructed. After that, 
all SPs (in Supplementary File 2) were discarded, and the 
new alignment of the data set without SPs, including 356 
informative characters of 492 total characters, was used to 
construct NJ, ML, and MP trees (Figure 2). Both the two 
data sets produced congruent phylogenetic topologies; 
however, for the high similarity of SPs among various 
LPs, the trees based on the data sets without SPs were 
more robust than those based on the data sets with SPs and 
thus were presented (Figure 2). Eight groups could be 
recognized and were designated according to the LP 



types in them (Figure 2). Among the eight groups, six 
groups, namely, LiP, MnP I, MnP II, VP, LiP and VP, and 
CII, were strongly supported (bootstrap values >90%) in 
all of the three tree-constructed methods; groups MnP III 
and MnP and VP well clustered together in all of the three 
methods (bootstrap values > 80%) and were with high 
bootstrap values > 90% in at least one method, so we 
treat them as groups. With regard to Phelebia radiata 
MnP and Trametes versicolor MnP02 PRL572, and 
Heterobasidion annosum pMnP5 and Pleurotus ostreatus 
pCIIl, they were strongly supported in the NJ method, but 
they did not cluster well in other two methods (bootstrap 
values < 50%o if cluster), so we treat them as separate LPs. 

Six groups contained only a single type of LPs and 
were found only in a single fungal order as follows: MnP I 
in Russulales, MnP II in Polyporales, LiP in Polyporales, 
VP in Polyporales, MnP III in Polyporales, and CII in 
Agaricales, while two groups had two types of LPs: MnP 
and VP found in Agaricales, and LiP and VP in 
Polyporales. It is notable that the LPs grouped neither 
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Phanerochaete chrysosporium LiPF RP78 
■Phanerochaete chrysosporium LiPC RP78 
Phanerochaete chrysosporium LiPJ RP78 
Phanerochaete chrysosporium LiPD RP78 

Bjerkandera adusta LiP 
Phlebia radiata LiPOl 
Phlebia radiata LiP03 
Phlebia radiata LiP04 
Trametes versicolor LiP07 
Trametes versicolor LiP 
Trametes versicolor LiPGIII 




VP (Polyporales) 



95/80/541 

100/100/100 L 

100/100/100 r Bjerkandera adusta VPOl ~| ^ 
'- Bjerkandera sp. VP J 

100/100/991 Heterobasidion annosum MnPOl B1256 

Heterobasidion annosum pMnP 
Heterobasidion annosum MnP02 B1169 
Heterobasidion annosum MnP03 B1257 
Trametes versicolor MnP 
Spongipellis sp. spMnPOl 
Phlebia radiata MnP 
Trametes versicolor MnP02 PRL572 
Trametes versicolor MnVi)5 KN9522 
Trametes versicolor MnPQ^ KN9522 
Trametes versicolor MnP04 PRL572 



LiP (Polyporales) 



MiiP 1 (Russulales) 



MnP II (Polyporales) 



—Pleurotus ostreatus MnPl 
- Pleurotus ostreatus MnP03 



100/100/100 p Pleurotus eryngii VP04 
100/100/99 Pleurotus ostreatus VP02 
\ 100/100/100 [ ~ Pleurotus ostreatus VPQl 
^ Pleurotus pulmonarius MnP05 

I Pleurotus ostreatus pMnP2 

[Pleurotus ostreatus VP04 
Pleurotus sapidus VPOl 
Pleurotus eryngii VP03 
Pleurotus eryngii VPOl 
Pleurotus eryngii VP02 

Trametes versicolor MnP^)\ KN9522 

Lentinula edodes MnPOl 

Heterobasidion annosum pMnP2 

Agaricus bisporus MnPOl 

Heterobasidion annosum pCII 

Trametes cervina CII 



MnP & VP (Agaricales) 



— Heterobasidion annosum pMnP5 
— Pleurotus ostreatus pCIIl 
-Pleurotus ostreatus pMnP3 
Trametes versicolor VP 

5Ganoderma fonnosanum LiP 
Ganoderma applanatum LiP 
luu, Ganoderma australe LiP _ 

-Laccaria bicolor pCll 

Dichomitus squalens MnPOl 

Dichomitus squalens MnP02 



LiP & VP (Polyporales) 



Phlebia radiata MnP02 

J Ceriporiopsis subvermispora MnPOl 

00/100/100 1 Ceriporiopsis subvermispora MnP03 

lOO/lOO/lOO i Ceriporiopsis subvermispora MnP02a 
' Ceriporiopsis subvermispora MnP02b 
[0 0/100/100 p Phanerochaete chrysosporium MnP03 
J ' — Phanerochaete chrysosporium MnP04 

' Ceriporiopsis subvermispora MnP04 

Phanerochaete chrysosporium MnPOl 
Phanerochaete sordida MnP03 
-Phanerochaete sordida MnPOl 
-Phanerochaete sordida MnP02 
-Phanerochaete chrysosporum MnP02 ATCC24725 
—Phanerochaete chrysosporum pMnP 
—Postia placenta CII 
- Coprinopsis cinerea CII 



MnP III (Polyporales) 



- Coprinellus disseminatus CII 



CII (Agaricales) 



-Malassezia globosa CcP 

— Ustilago maydis CcPOl 

— Coprinopsis cinerea CcPOl 

— Phanerochaete chrysosporiun 

Antrodia cinnamomea CcP 

Postia placenta CcP 



-Phanerochaete chrysosporium CII RP78 



Figure 2. Phylogenetic tree for 76 LPs and 6 CcPs of amino acid sequences without SPs. Topology is from NJ method, while bootstrap 
values (not less than 50%) from NJ, ML, and MP methods are indicated as percentages at the nodes. Groups LiP, VP, MnP I, MnP II, 
MnP and VP, LiP and VP, MnP III, and CII are delineated by vertical brackets at the right; meanwhile, the orders, which the LPs within 
each group belong to, are included in the brackets behind the group names. Six CcPs are set as outgroups. 
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entirely according to the LP types nor the fiingal phylo- 
geny of their original species. Even in the species from the 
same order, one type of LPs could have different origins. 
However, no LPs from the species of different orders 
grouped together 

Comparison of gene structures versus amino acid 
sequence alignments 

Gene structure analysis was performed in order to deeply 
probe the relationship among LPs. Certain minor varia- 
tions in intron positions were accepted, since that may be 
caused by indel or intron sliding (Stoltzfus et al. 1997). 
However, the predicted gene structure of Pc LiPH RP78 
differed from the previous study (Ritch and Gold 1992). It 
would be corrected if the last four nucleotides of the fourth 
exon and all four nucleotides of the fifth exon were 
transferred to introns, and the last eight nucleotides of 
the fifth intron were transferred to exons. However, the 
above adjustments would produce a noncanonical 3 '-spli- 
cing site. This abnormality may stem from errors in large- 
scale automated sequencing also found in other similar 
studies (Kawaguchi et al. 2007; Zhou et al. 2008), or it 
may be real as a similar noncanonical 3 '-splicing site has 
also been identified in other Ip introns. In either case, there 
is adequate information to suggest that the previously 
predicted gene structure of Pc_LiPH_RP78 was not cor- 
rect and it was consequently adjusted as described above 
and utilized in subsequent analyses. 

A total of 36 intron positions were identified and 
named from A to AJ (Figure 3). Several LPs were 
excluded from the comparison, because their gene struc- 
tures were unavailable (Table 1). This produced that only 
a single LP existed in groups VP and MnP II in the 
comparison, and thus it was unable to compare Ip gene 
structures in these two groups, whereas it would be pos- 
sible to do these in other groups, especially in groups LiP, 
MnP and VP, and MnP III with 17, 8 and 15 LPs, respec- 
tively, and thus obtain some general information. Whole Ip 
gene structures were not highly conserved and only one 
intron, K with a nearly identical phase 2, appeared in all 
groups. However, most introns coexisted in several 
groups, such as introns A in groups LiP, VP, MnP I, 
MnP II and MnP and VP, O in groups LiP, VP, MnP I, 
MnP II, MnP and VP, LiP and VP and CII, T in groups 
LiP, VP, MnP I and MnP and VP, and so on. Five introns 
appeared only in certain Ip genes: introns J in Tce CII, P 
in Led_MnP01, V and Z in Lb_pCII, and AF in 
Han_pMnP5. Remarkably, the introns within the same 
phylogenetic group, especially those simultaneously in 
the same species or even the same genus, were fairly 
conserved in both position and phase. For example, in 
group LiP, introns H, S, T, U, and AJ, intron R, and intron 
AG were exactly identical in each LP of P. chrysosporium, 
P. radiata, and T. versicolor, respectively; in groups MnP 



I, MnP and VP, and LiP and VP, each intron of LPs 
appeared in the same positions with identical phases; in 
group MnP III, each LPs possessed the same four introns, 
C, Q, Y, and AL 

Selective pressure acting on LPs 

The selective pressure on all LPs (see Supplementary File 
3 for their coding sequences) was evaluated using three 
pairs of site models. M8 indicated positive selection with 
CO = 2.38, and many positively selected sites were identi- 
fied by NEB, most of which were at the 99% level 
{P > 99%; Supplemental Table SI), but the LRT statistic 
of M7 vs. M8 rejected the hypothesis (Table 3). No 
unambiguous positive selection was thus found in the 
whole LP family. Subsequently, each of the eight LP 
groups was analyzed separately using the same three site 
models. Positive selection in groups MnP I, MnP III, and 
CII was ensured by MO vs. M3 (co = I.2I, P < 0.01) and 
M7 vs. M8 (« = 1.49, P < 0.025), MO vs. M3 {m = 1.56, 
P < 0.01), and MO vs. M3 (m = 19.89, P < O.OI), respec- 
tively, while in the other groups, positive selection was 
rejected either by an co value not more than 1, or by a P 
value not less than 0.05 in LRT statistics, following the 
rules of Yang (2000). 

Discussion 

In this study, we assembled the largest collection of LPs to 
date, explored their phylogeny, and examined the mechan- 
isms mainly caused the significant differences between 
white-rot and brown-rot fungi in degrading lignin. 

The most notable feature of the phylogenetic analysis 
was their topologic structures. The LPs of different types 
and/or from different species often clustered in the same 
group (MnP and VP and LiP and VP) with high bootstrap 
values, and vice versa (MnP I-III). This might be caused 
by the established standard for the phylogenetic topologic 
structure based on the information of overall amino acid 
sequences is different from the phylogeny of the LP types 
based solely on several key residues involved in substrate 
oxidation. It is the topologic structure rather than the LP 
types that reflects the bona fide phylogenetic relationship. 
Given the LPs in each group evolved independently and 
thus had their own unique ancestors, we could get some 
infomiation on the origin of LPs from the phylogenetic 
relationship: (1) the same types of LPs might originate 
from different ancestors independently, such as MnP from 
groups MnP I, MnP II, MnP and VP and MnP III, LiP 
from groups LiP and LiP and VP, and VP from groups VP 
and MnP and VP. (2) Even in the same species, different 
types of LPs might have not only one origin. For example, 
MnPs, LiPs and VP of Trametes versicolor (L.) Lloyd fell 
into groups MnP II, LiP and LiP & VP, respectively, but 
not assembled together So did LiPs and MnPs of P. 
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Figure 3. Comparison of the intron position and phase of Ips based on amino acid sequence alignment. Each intron is named from A to 
AJ according to its homologue, and its phase is represented by 0, 1, or 2. The LPs from the same groups were assembled together to 
visualize the comparison clearly, and their code names are the same as those in Table 1. 
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Table 3. LRT statistics and P value. 



Group 


MO vs. M3 (d.f. 


= 4) 


Mia vs. M2a (d.f. 


= 2) 


Ml vs. M8 (d.f. = 


2) 


2M 




P value 


"1 \ p 


P value 


O A f 


P value 


All sequences 








— 




-11044.875094 


— 


LiP 
















VP 


6.663864 




>0.10 


1.475886 


>0.10 


2.073648 


>0.10 


MnP I 


24.743772 




<0.01 


6.744184 


<0.05 


6.904412 


<0.05 


MnP II 
















MnP & VP 
















LiP & VP 
















MnP III 


220.201258 




<0.01 










CII 


56.385656 




<0.01 


1.777202 


>0.10 


2.001702 


>0.10 



chrysosporium, LiPs and MnP oi Phlebia radiata Fr., and 
LiP and VP of Bjerkandera adusta (Willd.) P. Karst. (3) 
Even the same type of LPs from one species might not 
cluster in the same group. Four MnPs from 
Heterobasidion annosum Fr. (Bref.) and two from T. ver- 
sicolor were scattered in the trees rather than clustered 
together within the corresponding groups MnP I and MnP 
II, respectively. This indicated that these MnPs either 
originated from different ancestors or experienced some 
relatively dramatic changes in amino acid sequences. (4) 
The LPs within the same groups only existed in an iden- 
tical order, although LPs from the same order could have 
multiple origins. This suggested that the origins of LPs 
were independent at the level of orders. In fact, during the 
screening of fimgal genome sequences, we found that nine 
taxa from other five orders outside Russulales, Agaricales, 
and Polyporales (Table 2) did not have any positive hits. It 
could be postulated that not all basidiomycetous orders 
had LPs, in despite that the genome of many basidiomy- 
cetous species are not sequenced at all. 

According to our searches of all available fungus and 
LP data, CII and MnP exist in the orders of Russulales, 
Agaricales and Polyporales, VP in Agaricales and 
Polyporales, and LiP only in the Polyporales (Table 1). 
The regularity of the distribution of LPs in the three orders 
was substantial, and thus it seems that with the aid of other 
evidence we could presume the evolutionary order of 
various LP types. For example, CII and MnP exist in all 
three orders and are more widespread than other types 
(Orth et al. 1993), suggesting that they were the earliest 
extant LP types. After them, VP, existing in two orders, 
appeared, and then LiP emerged in Polyporales. This 
evolutionary order met the need of white-rot fimgi to 
degrade the heterogeneous and recalcitrant lignin polymer. 
It could be postulated that in response to the appearance of 
lignin with high redox potential, the ancestor of LPs 
evolved from CcPs through the modifications of some 
key amino acid residues (Passardi, Bakalovic, et al. 
2007). CII without any obvious lignin-degrading activity 
and MnP degrading only phenolic lignin units were the 



primary groups of enzymes to evolve. Following them, 
VP and LiP with a capacity to degrade nonphenolic struc- 
tures, the most lignin units (Boeijan et al. 2003), emerged, 
and thereafter lignin could be degraded completely and 
efficiently. This conjecture is mostly based on the evolu- 
tionary order of the three basidiomycetous orders: from 
Russulales to Agaricales, and then to Polyporales. 
However, there is still ambiguity on the evolutionary 
order of the three basidiomycetous orders (Hibbett et al. 
2007), when various data sets and phylogenetic methods 
were used. So making sure the evolutionary order of the 
three orders would dramatically help us to elucidate the 
origins of various LP types. 

Besides the similarities of the amino acid sequences 
(directly reflected by the topologies of phylogenetic trees), 
Ip gene structures were also used to show the relationships 
among various LPs. The gene structures were, more or 
less, conserved within each group relative to those among 
different groups, suggesting that the LPs within the same 
groups did have their own common ancestors, despite that 
they had different enzyme types and were even from 
different species. The ancestors within each group either 
yielded a series of LPs with a single type or diversified 
into different LP types in various fiingal species. 

It seems easily to understand that the LPs from a 
common ancestor have an identical enzyme type, whereas 
more evidence is needed to explain how to produce not 
only one LP type from a common ancestor The topologies 
of groups MnP and VP and LiP and VP showed that the 
enzyme types could be transformed between MnP and VP, 
and LiP and VP, respectively. Furthermore, the transfor- 
mation has been successfully carried out in vitro by site- 
directed mutagenesis (Timofeevski et al. 1999; Mester and 
Tien 2001). Besides that, as previously mentioned, the LP 
lies on several residues involved in the process of oxidiza- 
tion. Therefore, selective pressures on LPs were analyzed 
to detect the dynamic of the transformations, namely, 
whether positive selection acted on LPs. 

The (o value and related parameters of LPs were 
estimated with PAML 4.2 using the CODEML program. 
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Positive selection was rejected in the whole LPs, but was 
doubtlessly ensured in the groups MnP I, MnP III, and CII 
under the used site models: MO vs. M3 and M7 vs. M8, 
MO vs. M3, and MO vs. M3, respectively. This indicated 
that the LPs in these three groups were with a potential 
bias to transfer to other LP types. The residues with the 
bias to mutation (Supplementary Table SI) were worth 
investigating further to show the functions of them in the 
oxidization reactions. Instead of tryptophan residue, tyro- 
sine residue (Y181) was found a catalytic residue in the 
LiP of Trametes cemna (Schwein.) Bres. (Miki, Calvino, 
et al. 201 1). Moreover, the chemically modified study also 
indicated that tyrosine could be considered as an effective 
catalytic site in the LiP of P. chrysosporium (Miki, 
Ichinose, et al. 2011). Therefore, positive selection also 
has a potential role in producing more unknown types of 
LPs. With regard to MnP type, only transferring to VP 
type was identified from the phylogenetic topology (the 
MnP and VP group). Through this transformation, LPs got 
a potential capacity to degrading high redox lignin struc- 
tures. However, it is hard to exclude the possibility that 
MnP type could transfer to LiP or CII types. Given that 
positive selection did not act on all LPs, but only on three 
of eight groups, adaptive evolution of LPs might only 
have occurred in a small fraction of fungal species that 
occupied a special niche. Large-scale ecological study will 
reveal why the positive selection occurred in these rather 
than in other species, and what was the driving force for 
the positive selection. Besides the above calculated evi- 
dence on the transformation of LP types, CIIs also occured 
in a few brown-rot fiangi, such as Ppl CII predicted in the 
P. placenta genome. Moreover, recently a functional LP 
closely related to VP was identified in a well-known 
brown-rot fungus Antrodia cinnamomea T.T. Chang and 
W.N. Chou (Huang et al. 2009). These phonomena indi- 
cate that the transformations of LP types could change the 
lignin degrading ability of wood-decaying fungi and even 
the rot types of them. In other words, if LP types in 
brown-rot fimgi were transformed to either MnP, LiP or 
VP from CII, the brown-rot fungi might also be changed 
to white-rot fiangi, and likewise if all LPs in white-rot 
fungi were transformed to CII type, the white-rot fungi 
might be brown-rot fiangi as well. 

Two other species of the order Agaricales, 
Moniliophthora pemiciosa (Stahel) Aime and Phillips- 
Mora and Schizophyllum commune Fr., lacked any type 
of LPs (Table 2). This might be considered normal for the 
former as it is a nonwood-decaying fiingus, but is surpris- 
ing in the latter case, as S. commune Fr. is a ubiquitous 
white-rot fungus with a worldwide distribution and its 
genome sequences detected in this study are quite com- 
plete (up to ~8.29x coverage). The presence of laccases, 
the second type of lignin-degrading enzymes, has been 
shown in S. commune (Hatamoto et al. 1999), suggesting 
that it contains a different system for lignin degradation 



and causes white-rot compared to the LP system in some 
other white-rot fiangi. Taking P. chysosporium with LPs 
and without laccases (Martinez et al. 2004) into considera- 
tion together, it seems that there are two distinct lignin- 
degrading systems, which convergently evolved in func- 
tion, and white-rot fungi obtained either one or both of 
them. Laccases also exist in insects, plants, and bacteria 
besides fungi, and they have multiple functions in differ- 
ent kinds of species (Nitta et al. 2002; Claus 2003; 
Arakane et al. 2005); however, recent evolutionary studies 
only focused on laccases from fiangi (Valderrama et al. 

2003) , and even Hoegger et al. (2006) including wider 
range of taxa did not elucidate the relationship among 
laccases from the four kinds of species. Therefore, phylo- 
genetic studies, reflecting to the origin and evolution of 
fungal laccases, were needed to explore whether the phy- 
logeny and lignin-degrading ability of fungal laccases 
entirely differed from those of LPs. 

The evolutionary pattern of multigene families has two 
models: concerted evolution and birth-and-death evolution 
as defined by Nei and Rooney (2005). As seen in the 
above discussion, sometimes LPs from the same species 
were not closer to each other than to those from other 
species in phylogeny, and there are at least twice wide- 
spread duplication events in their evolutionary history. 
One brought the ancestors of each group from the most 
primitive ancestor, and the other diversified the members 
of each group. Among the various LP types, besides MnP, 
LiP, and VP, there was also a potentially nonfunctional CII 
and perhaps also some other types derived from the above 
gene duplication events. Thus, the Ip gene family might be 
subjected to the birth-and-death evolution just as the 
Class-I of the plant peroxidases superfamily (Zamocky 

2004) . Similarly, in groups MnP and VP and MnP III, 
the LP clusters did not follow strictly the phylogeny of 
their original species. This indicated that the gene duplica- 
tion events probably occurred before the divergence of 
these fungal species or genera, and the duplicated LPs 
were passed down to newborn species or genera, respec- 
tively, after the divergence. The conservation of gene 
structures in each of the four groups also supports the 
theory that the LPs in these four groups were subjected 
to birth-and-death evolution. Conversely, groiap LiP con- 
tained 17 LiPs from four species, or say genera, belonging 
to the same order, Polyporales, and the LiPs from the same 
species/genus clustered together with strong supports. 
Furthermore, except for the introns shared in all genes of 
this group, there were also some species- or genus-specific 
introns. Thus the multiple LiP copies in one species or 
genus might be derived from an intraspecies/intragenus 
gene duplication event, which means that the LiPs in 
group LiP were subjected to concerted evolution. With 
regard to the groups MnP I, MnP II, VP, LiP and VP, 
and CII, the number of LPs in each one was too low to 
exactly infer the situation of gene duplication events. 
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much less to identify the evolutionary patterns in these 
groups. We can, however, conclude that the Ip gene family 
was subjected to the mixed process of both concerted and 
birth-and-death evolution. 

In summary, this study may provide essential evidence 
for elucidating the evolutionary dynamics of LPs and thus 
gaining insights into the alternation of white-rot and 
brown-rot fungi in their common lineage. LPs evolved in 
the mixed models of concerted and birth-and-death man- 
ners. They originated multi-times at the level of orders 
independently, and thus they did not exist in all basidio- 
mycetous orders, but only in several (three according to 
current study) orders. After the duplication events of 
initial LPs, which provided the raw materials for mutation 
events, various LP types might occur through mutations of 
several key residues driven by positive selection. This 
dynamics still acted on LPs from a small fraction of 
species, and might change LP types and even alter the 
rot types of wood-decaying fiingi. It is still hard to con- 
clude the evolutionary orders of various LP types. 
However, as the strict regularity of the distribution of LP 
types in Russulales, Agaricales, and Polyporales, it could 
be postulated as soon as we knew the evolutionary orders 
of the three basidiomycetous orders. These discoveries 
explain why white-rot and brown-rot fungi alternate in 
their lineage and will also provide a foundation for further 
research on the application of LPs in degrading persistent 
organic pollutants. 
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